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We derive a new first-order formulation for Einstein's equations which involves fewer unknowns 
than other first-order formulations that have been proposed. The new formulation is based on 
the 3-1-1 decomposition with arbitrary lapse and shift. In the reduction to first order form only 
8 particular combinations of the 18 first derivatives of the spatial metric are introduced. In the 
case of linearization about Minkowski space, the new formulation consists of symmetric hyperbolic 
system in 14 unknowns, namely the components of the extrinsic curvature perturbation and the 8 
new variables, from whose solution the metric perturbation can be computed by integration. 
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I. INTRODUCTION 

Many ways have been proposed to formulate Einstein's 
equations of general relativity in a manner suitable for 
numeric computation. In this paper we introduce a new 
first-order formulation for Einstein's equations. This sys- 
tem involves fewer unknowns than other first-order for- 
mulations that have been proposed and does not require 
any arbitrary parameters. In the simplest case of lin- 
earization around Minkowski space with constant lapse 
and vanishing shift, the system has the simple form: 
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Here the extrinsic curvature perturbation, a sym- 

metric tensor, and Xijk is a third-order tensor which is 
antisymmetric with respect to the first two indices and 
satisfies a cyclic identity, with the result that the system 
above is symmetric hyperbolic in 14 unknowns. 

Our approach applies as well to the full nonlinear ADM 
system with arbitrary lapse and shift. We work with the 
actual lapse, rather than a densitized version. In the 
nonlinear case the system involve 20 unknowns and can 
be written 

dohij = —2akij + 2hg(idj)¥ , 

-^dafimn = d[i{ak)^^n H • 

Here do := dt — b^di is the convective derivative, indices 
are raised and lowered using the spatial metric compo- 
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nents hij and the omitted terms are algebraic expressions 
in the hij, their spatial derivatives dihij, and the extrin- 
sic curvature components kij , and also involve the lapse a 
and shift components 6*. The fimm which depend on the 
first spatial derivatives of the spatial metric satisfy the 
same symmetries as in the linear case, and so represent 
8 unknowns. 

As is common, our derivation will start from the 
Arnowitt-Deser-Misner 3-1-1 decomposition ; see also 
Pj. The ADM approach introduces a system of coor- 
dinates t — xo, Xi, X2, X3, with t a timelike variable 
and the Xi spacelike for i = 1,2,3, and encodes the 4- 
metric of spacetime as a time-varying 3-metric on a 3- 
dimensional domain together with the lapse and shift, 
which are scalar-valued and 3-vector-valued functions of 
time and space, respectively. Specifically the coordinates 
of the 4-metric are given by 



500 



kb,h'^, 



5oi 



9i3 



Here a denotes the lapse, the bi are the components of the 
shift vector 6, and the hij the components of the spatial 
metric h. As usual Roman indices run from 1 to 3 and 
(h^^) denotes the matrix inverse to (hij). Let Di denote 
the covariant derivative operator associated to the spatial 
metric and set 



-Diibj), 



the extrinsic curvature. Then the ADM equations for a 
vacuum spacetime are 



dthij 



(2) 



dthj = a[Rij + {k\)kij - 2kiik\] + b^Dihj 



+ kuDjb^ + ki,D,b^ - D.Dja, (3) 
+ = 0, (4) 



D'k, 



D,k] = 0. 



(5) 



Here R denotes the spatial Ricci tensor, whose compo- 
nents are given by second order spatial partial differential 
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operators applied to the spatial metric components. Also 
indices are raised and traces taken with respect to the 
spatial metric. We have used the notation of parenthe- 
sized indices to denote the symmetric part of a tensor: 
D(ihj) '.= {Dibj + Djbi)/2. Later we will similarly use 
bracketed indices to denote the antisymmetric part and 
sometimes bars to separate indices, so, as an example, 

U[i\j\k]l '■= {Uijkl — Ukjil)/^. 

Equations l@J and (O , called the Hamiltonian and mo- 
mentum constraints, do not involve time differentiation. 
The first two equations are the evolution equations. A 
typical approach is to determine the lapse and shift in 
some way, find relevant initial data for h and k satisfying 
the constraint equations, and then to solve the evolution 
equations to determine the metric and cxtrinisic curva- 
ture for future times. The constraint equations may or 
may not be explicitly imposed during the evolution. For 
exact solutions of the evolution equation with initial data 
exactly satisfying the constraints, the constraints are au- 
tomatically satisfied for future times. 

The system of evolution equations for h and k is first 
order in time and second order in space. They are not 
hyperbolic in any usual sense, and their direct discretiza- 
tion seems difficult. Therefore, many authors |1,Q|EIS 
have considered reformulations into more standard first 
order hyperbolic systems. Typically these approaches in- 
volve introducing all the first spatial derivatives of the 
spatial metric components, or quantities closely related 
to them, as 18 additional unknowns. The resulting sys- 
tems involve many variables, sometimes 30 or more. In 
the formulation proposed here, we introduce only 8 par- 
ticular combinations of the first derivatives of the metric 
components. 

In the next section of the paper we present our ap- 
proach as applied to a linearization of the ADM system. 
This allows us to demonstrate the basic ideas with a min- 
imal of technical complications and to rigorously estab- 
lish the relationship between the new formulation and 
the ADM system. In the third and final section of the 
paper we carry out the derivation in the case of the full 
nonlinear ADM system. 



II. SYMMETRIC FORMULATION FOR THE 
LINEARIZED SYSTEM 

We linearize the ADM equations about the trivial so- 
lution obtained by representing Minkowski spacetime in 
Cartesian coordinates: hij = 5ij, kij — 0, a = 1, 
bi — 0. Consider a perturbation given by hij — Sij + jij, 
kij = Kij , a = 1 + a, bi = Pi, with the jij, Kij, a, and f3i 
supposed to be small. Substituting these expressions into 
the ADM system and ignoring terms which are at least 
quadratic with respect to the 7^ , Kij, a, and Pi and their 



derivatives, we obtain the linear system we shall study: 

dtjtj = -2k,j + 2d(^,Pj-), (6) 

dtKij = {Pj)ij - didja, (7) 

(^7) ■ = 0, (8) 

gj^^^ - d,K^ = 0. (9) 

Here P7 is the linearized Ricci tensor, with components 
given by 

(^^7)^ = ld^d'-/i, + ld,d'ju-^d'dn,,-^d,d,jl. (10) 

In these expressions, and in general when we deal with 
the linearized formulation, indices are raised and lowered 
with respect to the Euclidean metric in M.^, so, for exam- 
ple, 9* and di are identical. 

In order to reduce the linearized ADM system to first- 
order symmetric hyperbolic form, we first develop an 
identity for P7 valid for any symmetric matrix field 7. 
It is useful to introduce the notations § for the the 6- 
dimensional space of symmetric matrices and T for the 
8-dimensional space of triply-indexed arrays (wijk) which 
are skew symmetric in the first two indices and satisfy the 
cyclic property Wijk + Wjki + Wkij = 0. 

We define the operator M : C°°(R^ §) C°°(M3, K^) 
by {Mu)i = d^Uii — diu\. (Note that the linearized mo- 
mentum constraint, is simply Mk = 0.) From the 
definition H10|l . we have 

(P7)y=-9'9[/7,]. + ^a,(Af7),. (11) 
Now, for any vector- valued function v, we have 

Applying this identity to (|ll|l with v — M7 we get 

(P7)., = -^\^^^J]^ + (M7)[;'5,],] + \d\Nh)A,. (12) 

Define operators L : C°°(M3,§) C°°(M3,T), and 
L* : C°°(K3,T) ^ C°°{R^,S) by 

{Lu)iji = d[iUj]„ {L*w)ij = -d^wi(.ijy 

One easily verifies that operators L and L* are formal 
adjoints to each other with respect to the scalar prod- 
ucts {u, v) — J UpqvP'' dx and {z, w) = J ZpqrW^''^ dx in 
the spaces C°°{R^,E>) and C°°(M3^t) respectively. Intro- 
ducing 

>^l3^ = + iMl)llSJ]^] (13) 

we can then restate (|12ll as 

(P7)„- = V2d'Xij, + ia'(M7)i<5„-. 
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(The reason for the factor of \f2 wih become apparent 
shortly.) Taking symmetric parts, the last equation be- 
comes 

{P^\, = -V2{L*\\, + ia'(M7)z%. (14) 

We are now ready to introduce our first-order sym- 
metric hyperbolic formulation. The unknowns will be 
K e C°°(IR3,S) and A e C°°(R3,T), so the system has 
14 independent variables in all. Substituting 114() in ((7|) 
and noting that d^{M^)i = {P"i)\ = by the linearized 
Hamiltonian constraint © , we obtain an evolution equa- 
tion for k: 

dtKij = -V2{L*X)ij - didja. (15) 

To obtain an evolution equation for A, we differentiate 
((T^ with respect to time and substitute © to elimi- 
nate 7. Simplifying and using the linearized momentum 
constraint © we obtain 

dtXiji = V2{LK)iji - Tiji, (16) 

where t can be determined from '■— d{iPj) by 

Ti,, = \/2[(ie/?)ij, + (Me/3)[,(5j],] 

= ^(9»9[,/?j] + d^dy^di^k, - 9"9[„/?,]<S,0- 

Equations H15(l and (|16|l constitute a first-order symmet- 
ric hyperbolic system (this is clear, since L and L* are 
formal adjoints). It follows that (see, e.g., j2j § 7.3.2), if 
the laspe and the shift are sufficiently smooth then for ar- 
bitrary initial values k(0) and A(0) belonging to i?^(R'^), 
there exists a unique solution to system H15|) and 1)16(1 
with components in iJ^((0, T) x M^). 

The Cauchy problem for the original linearized ADM 
system consists of the equations ©-© together with 
specific initial values 7(0) and k(0). The foregoing deriva- 
tion shows that if 7 and k satisfy the ADM system and 
A is defined by (|13|l . then k and A satisfy the symmetric 
hyperbolic system p5|) . (|16|l . Conversely, to recover the 
solution to the ADM system from H15|l . ((16|l . the same 
initial condition should be imposed on k and A should be 
taken initially to be 

Ay»(0) = -^{L-^m^^ + (M7(0))[„5,],]. (17) 

Once K and A are determined, the metric perturbation 7 
is given by 




as follows from ijHJi. 

Theorem 1. Let the lapse perturbation a and shift per- 
turbation (3 be given. Suppose that initial data 7(0) and 
k{0) are specified satisfying the constraint equations ijSJl, 



(|^ at time t = 0. Define A(0) by (|17() . and determine k 
and A from the first- order symmetric hyperbolic system 
(Cni, Finally, define 7 by (dJ). Then the ADM 

system ©-(PJ is satisfied. 

Proof. Equation © follows from H18I) by differentiation. 

Next we verify the momentum constraint To do 
so we will show that /i := Mk satisfies a second order 
wave equation, and that fi{0) = 9t/i(0) = 0. Indeed 
/i(0) = by assumptions. To see that dtn{0) = 0, we 
apply the operator M to (|15|) and use the fact that M 
annihilates the Hessian didja for any function a. There- 
fore dtfJ. — —^/2ML*X. Using (|14|l and the assumption 
that 7 satisfies the Hamiltonian constraint at the initial 
time, we find that 

^t^l^{o) = (MP7(o)), - -^d,d'{M^{o))i - 0. 

To obtain a second order equation for /i, first we differ- 
entiate ((15() in time and substitute p6|l to get a second- 
order equation for n: 

d^Kij = —2{L* LKij) — didja. 

Here we have used the fact that L*t = 0. Apply M 
to the last equation. Using the identity {ML*Lk) — 
— 9'9(i(AfK);)/2 and the fact that M annihilates Hes- 
sians, we find that /i := Mn satisfies the second-order 
hyperbolic equation 

This is simply an elastic wave equation. Since /i(0) = 
dtfi{0) = 0, fi vanishes for all time, i.e., the momentum 
constraint is satisfied. 

Now {PK)i = d\MK)i = 0. Moreover, P apphed to 
6/9 is identically zero. Therefore if we apply P to l|18|) 
and take the trace, we find that {P^)\ — (^7(0))^, which 
vanishes by assumption. This verifies the Hamiltonian 
constraint ©. 

It remains to verify Q which, in view of H15|) comes 
down to showing that \f2L*\ = — P7. Since we have 
verified the Hamiltonian constraint, this will follow if we 
can establish (|14|l , which is itself a consequence of H13|l . 
We used (|13|) at time t = Oto initialize A, so it is sufficient 
to show that 

^tXlj^ = --^dt[{L'y)ij^ + (A/7)[,(5j]i]. 

This follows directly from lfTC|) and □ 

We conclude this section by computing the plane wave 
solutions to the hyperbolic system That is, we seek 
solutions of the form 

Kij = Kijf{st - n^Xa), Xljt = Xijifist - n^Xa), 

where the real number s gives the wave speed and the 
unit vector Ui the wave direction, the polarizations kij 
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and Xiji are constant, and the profile f{t) is an arbitrary 
difFerentiable function. 

Substituting these expressions into (0), we get 

Skijf'{st - n^Xa) = -\/2n'A;(jj)/'(si - n^-Xa), 

sXijif'{st - n'^Xa) = -V2n[iKj]J' {st - n'^Xa), 

and so the system reduces to the following linear eigen- 
value problem: 

skij = -\/2n'A;(jj), (19) 
sXiji = -V2n[iKj]^, (20) 

with the wave speed s as eigenvalues and the pairs 
(kij,Xiji) as eigenvectors. The eigenvalues of this sys- 
tem are (multiplicity 4), ±1 (each multiplicity 3), and 
±l/\/2 (each multiplicity 2). To verify this and describe 
the eigenvectors we introduce a unit vector rrii be a per- 
pendicular to Hi, and set li = ef'^nafrib to complete an 
orthonormal frame. Then the following solution to the 
eigenvalue problem can be checked by direct substitution 
into 

s = 0: (0, mff^jjmi), {O^myilj^^k), (nin-,,0), 

(0, 2m[ilj]ni + n[ilj]mi - n^imj^U); 

s = ±1: {klj - m^mj ,^V2{n[ilj]k - n^inij^mi)), 

{kruj + milj,^V2{n[ilj]mi + n^imj^U)), 
{klj + niimj ,^V2{n[ilj]k + n^inijimi)); 

s^±l/V2: {n(ilj),^n[ilj]ni), 

{n(imj-j,Tn[imj]ni). 

III. DECOMPOSITION OF THE FULL ADM 
SYSTEM 

In this section we develop a first order formulation of 
the full nonlinear ADM system analogous to that devel- 
oped for the linearized system in Section ^ We con- 
tinue to assume that the underlying manifold is topolog- 
ically and view the ADM system as equations for the 
evolution of a Riemannian 3-metric /i on R^. Thus hij, 
1 ^ *jj ^ 3 are the components of the spatial metric. 
They form a positive-definite symmetric matrix defined 
at each point of R'^ and varying in time. Indices on other 
fields are lowered and raised using hij and the inverse 
matrix field h^^ . 

For the components of the Riemann tensor we have 

Rijkl — dj^ikl jkl ^" h jkmX^ iln ^ ikmT^ jItl) 

where the Cristoffel symbols are defined by r^T; = 
(dihjk + djhki — dkhij)/2. The components of the Ricci 
tensor are given by Rij — hP'^Rpigj, which yields 

= ^h-P'^idpdjhiq + didphqj - dpdqh^j - didjhpq) 

^ lb II ips-^ qjr pqs^ ijtJ- 



We define a second-order linear partial differential op- 
erator P : C°°(R3,R3><3) ^ C°°(R3,R3x3) by 

{P^h = l[dp{h'''d,Uq,) + d,{hPWpUjq) 

- dpihPldqU,,) ~ d.ihPWjUpq)]. 

Note that in the case hij = Sij , P coincides with the lin- 
earized Ricci operator introduced in Hl()|l . In the nonlin- 
ear case, P is related to the Ricci tensor by the equation 

i?.y-(P/i).,+c3, (21) 

where 

Clj -^[{dphP'^){d,h,q) + id,hP''){dphq,) 

~ {dphP'i){dqK,) - {d,hP'i){d,hpq)] 
-r It 11' y-L ips^ qjr pqs^ tjr J: 

is a first-order differential operator in hij. 

Define an operator M : C°°(R3,S) ^ C°°(R3,R3) by 

(MU), = 2hP''d[pU,]q. 

This formula extends the definition of the linearized mo- 
mentum constraint operator introduced in the previous 
section. Up to lower order terms the momentum con- 
straint (|SJ) is given by the vanishing of Mk: 

DPhp~D,kP = 2hP''D[ph]q = (Mfc),-2/iP«r;[pfc,], (22) 

(where T^p := h^'Tqps). 

Finally, we introduce operators L : C°°(R^,§) — > 
C°°(R3,T) and L* : C°°(R3,T) ^ C°°(R3,§) by 

{Lu)iji = d[iUj]i, iL*v)ij = -hqihrjdpvP'^'^''\ 

As in linear case, the operators L and L* are formal 
adjoints to each other with respect to the scalar prod- 
ucts {u, w) = J UijW^^ dx and (w, z) = J Vijiz^^^ dx on the 
spaces C°°(R3,S) and C°°(R3,T) respectively. 
Finally, we introduce new variables 

flmn = ——J=[[Lh)lmn + (Af/l) , (23) 

and develop the analogue of the identity (|14|l . 

Lemma 2. The following identity is valid for the Ricci 
tensor 

R,, = -V2{L*f),, + ^dp[hP'^{Mh)qh,,] + 4\ (24) 
where 

4 = 4j) + hq^hrJ[^p{h'''^hn]h''\{Lh)nran) 

is first order in hij . 
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Proof. The formula H24|l is a consequence of the identity Here we used the fact that the Lie derivative tfDikij + 



where 



.3) - 



To prove the identity (|25|l we note, first, that operator P 
can be rewritten in terms of operators L and M as 

(Pu)y = ~dp[hP'i{Lu)qj^] + \d,{Mu), 
which yields 



kiiDjb'' + kijDiU = b^dikij + kudjU + kijdib'. We treat 
the second term on the right-hand side of H24|l using the 
Hamiltonian constraint Q). Now, 

i?- = ^h''hP'^{dpdjhiq + didphqj ~ dpdqhij - didjhpq) 

-Til II IL yi^ ips^ qjr ^ pqs^ 2jr J 
= h'^hP'i{d,dphq,-d,d,hpq) 

+ h^ihPih^^(r r . _r r - ) 

-r 11/ IL y± ips^ qjr ^ pqs^ ijr J 
= d,[h''W"'{dphq,~d,hpq)] 

-[d,{0hP'i)]{dphqj-d,hpq) 

+ hPi h^'^ (T r -r r- ) 

-r 1^ ips^ qjr ^ pqs^ IJT } 

where 



= -dp[W"^{Lu)q(,j)\ + -d(,{Mu)j), 
and, second, that according to the definition of L* , 

{L*V),J = -hq.KjdpihP'h'i^^h^^Vli^^^)) 
- -dpih^'^Vq^,^)) 

- hq,hrj[dp{h1"'h'-"Wvn,nn)- (26) 

To derive the identity H25|l we substitute —y/2fqij for Vqij 
in (|26() . The first term on the right-hand side of (|26|l then 
becomes 

-dpihP''[iLh)q^,,) + l{Mh)qh,j - ^(M/l)(./l,),]) 

= -dp[hP'^{Lh)q^,^)] ~ ^dp[hP'^{Mh)qh,,] 

+ ldp[hP'^iMh)^,h,)q] 

= -dp[hP^Lh)q^,^j] - ^dp[hP''{Mh)qh,,] 
= iPh\,,^~^dp[hP^Mh)qh.,]. 

The substitution of — \/2/imn for i^imn into the second 
term of the equation (|26|l gives the term c^j precisely. 
The rest of the proof follows from the identity (|21|l . □ 

We now proceed to the derivation of the new formula- 
tion of the ADM system. In Q we substitute do + Udi 
for dt and replace Rij with the right side of the (|24|) to 
get 

5o% = -V2a{L*f),, + ^adp[hP''{Mh)qh,,] + c%, (27) 
where 



S = -{d,(h'^hP'i)\{dphq, - d.hpq) 

~ \^ ips^ qjr pqs^ ijr }• 

Hence, due to the Hamiltonian constraint, 

dp[hP'J{Mh)qhrj\ 

= {dp[hP'>{Mh)q]}h,j + hP"{Mh)q{dpKj) 

= [kpqkP'i - (fcP)2 - C^)]/ly -f hP%Mh)q{dph,,). 

Combining all lower order terms into an expression Bij^ 
first order in we reduce (|27|) to 



-V2a{L*f),,+B,j. 



(28) 



This is the first evolution equation of our system. 

The second evolution equation will be obtained by ap- 
plying ^0 to the definition of / H23|l : 

dofhnn = --^{do{Lh)ljnn + do[{M h)[ih^]„]} . 

First, we note that 

do{Lh)lran = (Ldoh) 

+ Imbmdsh^n) - {d^bWshin)]. 
Using the fact that {Mu)i = 2hP'^{Lu)piq, we then get 
do{Mh)i = {Mdoh)i+2{dohP'')iLh)piq 

+ hP1[{dpb'){d,hlq) - {dlb'){dshpq)]. 

If we use this formula to compute do[{Mh)ih„m] and then 
antisymmetrize in I and to, we obtain 

doflmn = - -^[iLdoh)l,nn + {Mdoh)[ihm]n] 



V2(ao/i^«)(L/i)p[,|,|V]„ 



(29) 
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where 

- ■^hP''[{dpb''){dshq[i)h„,],, ~ {dshpg){d[ib')h„,]^]. 

Next we use Q to relate the terms in (|29f) involving 
doh to the extrinsic curvature k. For the Lie derivative 
of the metric, we have 

2D^,bj-, = b^dshij + 2h,(,dj)¥. 

Using this l(2Jl becomes 

d^hij — —2akij + 2wij, (30) 

where Wij :— h^j^^dj^b'^. Using the Leibniz rule we can 
then verify that 

Substituting these expressions in (|29|l we obtain 

dofimn = \/2[L(aA;)]i„i„ + \/2[M(afc)][i/i„]„ + c[^j„, (31) 
where 

Cllnn = cf^„- V2[(Lw);™„ + (Afw)[,/l„]„] 

-2y2[aF«-u;f«](L/i),[,|,|/i„]„ 

+ ^/2{Mh)[iakm\n - V2(M/l)[;W„]„. 



The final step is to invoke the momentum constraint to 
simplify the second term on the right-hand side of l|31|l . 
Indeed, since the right-hand side of H22|) vanishes, 

[M{ak)]i = aiMk)i+2hPiid^a)ki^g 

= 2a/i^«T^[pfc;],+2/iP«(5[,a)fci],. 

Substituting this in H31|) we obtain the desired second 
evolution equation: 

V2[L{ak)] + Ci (32) 

where 

-I- hP'^{d[pa)ki]ghmn - /i'"'((9[pa)fc,„],/i/„]. 

The two equations l|28() and (|32|l constitute a first order 
system for the unknowns fcy and fimn- This system is 
coupled to the ordinary differential equation H3U|) through 
the terms Bij and Cimn which are algebraic combinations 
of hij , di hij , kij , the lapse a and the shift b and their 
spatial derivatives. The foregoing derivation shows that 
if h and k satisfy the ADM system then h, k, 

and / satisfy the system (jSOl, (I2H1), 
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